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OPERATOR 

Abstract 

Some consequences which follow from the periodicity assumption for spectral density of Wilson-Dirac 
operator are studied. Such an assumption allows to obtain simple representations for quark propagator, 
which reveals an important role of m <-> —m symmetry. It is argued that this symmetry is restored when 
the mirror fermion mass m r tended to infinity. The constrains on zero modes of Wilson-Dirac operator 
in a toy model approximation are also discussed. 
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1 Introduction 

The spectrum of Dirac operator is closely related to a number of important aspects in finite 
temperature QCD. The behavior of spectral density p(X) around A ~ for massless Dirac operator 
is of particular physical interest being responsible for breaking the chiral symmetry [[l]. The low 
lying eigenvalues of staggered fermion operator are also used to extract the topological content of 
background gauge field configurations (|]. 

Despite steady progress, the numerical simulation results for p(A) remain inconclusive. As it is 
known from general theoretical considerations (see e.g. || ||, ||), the eigenvalues of the normalized 
Dirac operator are located in a compact area. However, as it pointed out in || perturbative 
calculations show even an increase of state density at large A. Indeed, for A 3> Aqcd the spectral 
density is regarded to be insensitive to gluon vacuum fluctuations and behaves in the same way 
as for free fermions: 

This free term introduces quadratically divergent piece in (iptp). However, to get a truly non- 
perturbative quark condensate which is the order parameter of the spontaneous chiral symmetry 
breaking, this perturbative divergent part should be subtracted. As a result, the quark condensate 
is related to the region of small A: A ~ m -C Aqcd ||- 

Some numerical results j?], |^, P may be taken as an evidence in favor of compactness of Dirac 
operator spectrum (see Fig. 2,4,5 and 6 of ref. M). Indeed, the spectral density disappears outside 
certain region T\ |u^], that for a particular case may be specified as: |A| < 2Aqcd (see Fig. 2 



of ref. H), or for |A| < 3/C (see Fig. 2 of ref. 1 1 1 ). The domain wall fermion approach leads to 



a similar picture (see Q Fig. 12). In case of (l+l)-dimension it was proved [l4|, [15], that 
the spectrum was contained inside the compact area (circle of radius 2 in the complex plane). 
Moreover, compactness more often than not may be treated as implicit periodicity, so it is hard 
to avoid the assumption that the spectral density may appear a periodic function of A. Indeed, 



the lattice regularization for Dirac operator leads to 

3 

D (A) = -f(dv - iA v ) -rn^J2 D • ( 2 ) 

where 

D W*'* = r —^- ( x ) 6 *,*'-» + r -^ JLU t (*') Sx^+u - (r + Slma) 8 X , <X (3) 

and r = tanh a T m r Wilson parameter. Let us assume that the fermion part of the action — Sp = 
lndet D can be expressed in terms of Polyakov loops £1 (x) = Yit^a 1 ^° x ) ' as ^ can be done 
in numerous particular cases for the gluon part of the action, and det D — f ((f) . In this case one 
may expect that the Dirac matrix detl? (ip) (as well as (tpip)) became periodic in <p 

det D (A ) = det D (A a + 2ttT) , (4) 

with 

(p = argQ = a T N T Ao = ^-. (5) 

Therefore, one may suspect that this condition, together with the periodicity along t - axis, may 
induce the periodicity of D in the complex plane of m and, hence, in eigenvalues A. In other words 
it would not be too surprising if 

A„ (<p + 2vr) = A„ (<p) + 2vr6; 6 = const. (6) 

Indeed, if to take the free Dirac operator as an example | fli"| , the spectrum in the phase 
argft = is given by: A 2 = k 2 + ((2n + 1)ttT) 2 and A 2 = k 2 + ((2n + 1/3)ttT) 2 for argfi = 2vr/3. 

More rigorous and convincing evidence in favor of the periodicity of Dirac operator spectrum 
in general case was suggested in |n| (see also [[l7| ^). Indeed, Dirac operator can be presented 
(§ as 

aD = 1 - l5 e(H)), (7) 
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where e(H) some Hermitian matrix and H = 75ft — 7 M (<9 M + iA^)]. Therefore, if D satisfies the 
Ginsparg- Wilson relation |20| 

D^o + I5D = aDj 5 D, (8) 

then the operator j 5 e(H) = 1 — aD is unitary 

(1 - aD ) f (1 - aD ) = (1 - a lb D l5 ) (1 - aD ) = 1. (9) 

This implies that the spectrum of D lies in the complex plane [H], [l8| on the circle, (1 — e l9 )/a 
with 9 6 [— 7r,7r]. Considering the above, a Banks-Casher-like relation Q will be 

(^) = (i?- 1 (a, a)) ~ a f P (g) w d6. (10) 
u /_ 7r 1 + am — e w 

The periodicity of Dirac spectrum p (6) in 6 may be, as a matter of course, treated as a lattice 
artefact Indeed, assuming that the major part of eigenstates is concentrated in the area 9 2 < a 2 , 
||l8[ after redefinition of the integration variable 9 — > sa, and lifting lattice regularization (a — > 0) 
one comes to the standard Banks-Casher relation 

(to) = (D-* ( X , x ))~r a r^^ds * r fHds. (ii) 

Though, one should also consider the possibility that e.g., A 2 < 1/T 2 , with A = 9T. In this 
case the spectrum remains invariant under the shift A — > A ± 2nT. The Banks-Casher relation, 
of course, remains valid in all cases, however, premature lifting the lattice regularization in ( |ll| ) , 
leaves the conjectural periodicity hidden. 

Certainly, the arguments given above, will not suffice as such and we must admit, that the 
assumption (^|) looks much more controversial then the evidences given in favor of compactness. 
Below we briefly sketch the results [^IJ obtained in the approximation, where all terms proportional 
to l/£ have been discarded]^] in a Hamiltonian limit a T — > (£ — > 00) . The spectral density 
computed in such an approximation is actually periodic. This reassures us a bit and allows us to 

1 The anisotropy parameter £ is defined as £ ~ a a /a r , where a T and a CT are temporal and spatial lattice spacings, 
respectively. 



assume that the periodicity may appear to be a general feature of the spectral density of Dirac 
operator. Hereafter we consider some consequences which follow from such an assumption for a 
more general case. 

2 Toy model approximation 

In a toy model approximation the fermion determinant can be computed analytically and presented 



in a closed form. We choose the form suggested in |22f] (see also 23 ) for the fermion part of action 
-$ = n f L 3 jy x ,D° x , x i, x + Xjk^'xTM , (12) 

\ x ^ n— 1 x / 

where rif refers to the number of flavors^] and D x , x are given by (||). 

Inasmuch as in the suggested approximation we omit the terms of l/£ order, the 'toy' action 
will be simply a?'Yl l jjj x D® xx ,^ x r. The projectors 1= ^ 7 " divide the bispinors if} into two components 
^(±) _ i±7o ^ each including only one two-component spinor. Therefore, taking into account 
that (S x -fi,,x') = S x ,x-p,> = S x +p,,x', and presenting 



r ± 7o l + 7 r±l r t 1 1 - 7o 



(13) 



2 2 2 2 2' 

we may rewrite ( |l2|) as 

-S F = v/T^ + V^A^vif 5 ) > (14) 

where 

A XK / = d xx < I <ty,t-i ^ ~ ~^ g ^t'.t I • (I 5 ) 

By gauge transformation all Uq (x, t) matrices may be diagonalized simultaneously: Uq (x, t) 
5 au UQ (x, i) , therefore A xx i in jl^ ) is simply a set of TV matrices AT T x iV T . Hopefully, the com- 
putation of det A, in the considered extremely simple case can be done straightforwardly and after 
the integration over ipi^ we get 

N 

—S e J^ = rif ^ In det A^Aq, = n/ln J^J (z r — cos^q) — (— l)^ cos<^ Q ^ , (16) 



We consider degenerate quarks. 



with 

2=(-l) cosh—; z r = (-1) cosh - , (17) 

where B = for periodic and _B = 1 for antiperiodic border conditions on fermion fields in a 
temporal direction. Thereafter we shall consider the standard case B = 1 and even values N T , 
because the generalization on B = and odd N T is quite evident. 
It easy to see from ( |l6| ) and ([l7]) that ( iptp) may be presented as 

(M) = 2n f sinh— < — i \ (18) 

v 7 J T \ cosh ^ + cos tp I K ' 

. m + 2m r / 1 \ 

+ In t sinn ( ^ ) , 

J T \ cosh m±phL + cos y>/' 

where m r is the mirror fermion mass connected to Wilson parameter by the relation a T m r = 

arctanh r. 

3 General case 

Recall that the fermion Green's function in the external gauge field is given by 

lit L /\ n 

n 11 

where \ n {x) eigenvalues of the massless Dirac operator and u n are corresponding eigenfunctions. 
The spectrum of the massless Dirac operator is discrete and enjoys the chiral symmetry. In 
particular, for any eigenfunction u n {x), the function 

U n = l°U n (20) 

is the eigenfunction as well, with the eigenvalue 

K = -A„. (21) 

Therefore 

p(X)=p(-X) (22) 



and eigenfunctions occur in pairs with opposite eigenvalues, except for zero modes^J . In the limit 
V — > oo, the level spectrum becomes dense and we can trade the sum in (|l9|) for the integral and 
get 

where 

In case when spectral density is compact and periodic with the some period 

p(X) =p{\ + 2kirO) (25) 



we can write 



/-27r(fe+l) (\\ 

m m = m E / ^% dX (26) 



or 



m m = m e / i , , „ u ^ dX - w 



k— — OG 

Presenting (pi 5.1.25(3)) 



o m 2 + (A + 2Trkey 



™ 1 sinh § (28) 

fe ±- oo (A+2fc 7 r) 2 + (f) 2 cosh f- cos A' 



we get 



<*>-* S-fsaffe* (29) 

So we come to the conclusion that ( |l8| ) is nothing but a particular case of ( p9| ) for the periodic 
spectral density p (A) with the period O = T. It easy to check that it leads to the expression for 
(^V) that coincides with ([l8]) obtained in a toy model approximation. 



3 Generally, symmetry 1221) does not mean chiral invariance. Indeed, the eigen functions of Dirac operator are 
the same for m = and 0. Therefore, the symmetry (122]) is preserved for arbitrary m, though only massless 



Dirac operator is invariant under transformation (pc|). 



If one claims 



p(X) = p(Ai; A2A3) = p(Ai + 2 7 rfc 1 T; A 2 + 2irk 2 T; A 3 + 2irk 2 T) (30) 



and 



^ = 27rfc T, (31) 

where all kj are integer, one can easily get an expression for (f/'V') which coincides with the result 
obtained in for SU (3) gauge group case^J . So we come to the conclusion that the 'toy' model 
gives the results that may appear reasonable in the general case (£ ~ 1) if periodicity assumption 
happens to be allowable. 

4 Some formal consequences of discrete symmetries 



A very formal expression of the symmetry (|20j) , (|22j) can be found if we take into account the fol- 
lowing. As it is known the transformation under CPT reversal can be presented as a combination 
of the unitary transformation^] ip (t, x) — » j^tp (t, x) and inversion ip (i, x) — > ip (— t, — x). Under 
unitary transformation as well as under inversion the massless Dirac operator 2) only changes the 
sign. Therefore, the result of the unitary transformation of J) — to can be mimicked by the formal 
substitution m — ► —to. It is easy to check that p (A) «-> p (—A) invariance also leads to to <-> — m 
symmetry 

E (to) = (^V) = r ^KdX = [°° ^^-d\ = — E (-to) . (32) 

To be specific, we confine ourselves to the case of SU (2) gauge group. Taking into account 
that 

— — i ^ = 2j>l) 2j x, feOe-W +1 "", (33) 

cosh a; ± cos 5 ^— ' 

2 j=0 



4 Of course, specific form of periodicity in ( [?o| ) was chosen intentionally to achieve such coincidence. 

5 This invariance under chiral rotation ip — > exp {2756*} i/> rotated by an angle 9 = it is, in fact, the remnant of 
chiral symmetry of the massless theory. 



where \j if) — 2 are the characters of irreducible representations of SU (2) group, we 

can obtain from 



= 4n / sinh^e-( 2 -'- +1 ^l (-1)" ( Xj ) + (34) 

3=0 

4n /S inh^±^f e-lW\*?*\ {-!)« , 

3=0 



so the breakdown of the symmetry (|22|) (and consequently (|32|)) by the second term in ( p4[ ) for any 
m r ^ became apparent. Moreover, Sp posses an undesirable invariance under the interchange 

m <-> — to — 2m r . (35) 

This invariance is inherent neither for Dirac operator nor for the spectrum ( p3| ) . It the limit a T — > 
such symmetry is cognate to the invariance under the hopping parameter interchange: K —f —K, 
which plays an important role in the emergence of conjectural Aoki phases, widely discussed in 
literature (see e.g., p5, 26, p7|). 



Although the mirror fermion input does not disappear in a T — > limit, almost all mirror terms 
exponentially vanish for ^5f — > oo. The only surviving term j = 0, however, introduces finite 
contributions into Otpi(}) which breaks m — * — m symmetry. Let us consider a very special case 
m r — ► oo , in other words, put exactly r — 1 in the action (^) /rom t/ie beginning. 

It is easy to see that in this case 

A x ,, x ^ x , x (U (x, t) 5 t , t /_ 1 - (ma T + 1) £ t , t ,) ~ S x , x (U (x, i) 5 t , t /_ 1 - e a - m 8 t , v ) , (36) 



therefore 



detA x> = JJ^+e^), (37) 



and 



exp {— Sj?} = n/ det A det A* = nf Y[ (cosh — + cos</? a ^ (38) 
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where e i(Pa are eigenvalues of the Polyakov matrix Q, a p = & a pe %Va . In this case (at least in zero 
order in l/£ ) mirror fermion term does not appear at all and terms leading in l/£ are invariant 
under the change m — > — m . The mirror fermion contributions in (pTi^are shifted into the area of 
infinite m r and undesirable symmetry ([35]) is expired. 

It is interesting to note that, after 'harmless' substitution 1 + a T m ~ e a ^ m the periodicity in 
A may appear in a general case. Indeed, we may expect that in the limit (a T — > 0; N T — > oo) the 
mass variable will survive only in combination [e a " r7n ) r9 = e^ q , where q is some dimensionless 
number. Thereby we get periodicity m — > m + 2iriT/q . Taking into account that eigenvalues A of 
the operator X) = 7^ [d v — iA v ) defined by det (D — iX) = are related to eigenvalues A m of the 
operator 2) — m by A m = A + im we also get (as can be seen e.g., from (p3[)) periodicity in A 

A -> A + 27rT/q. (39) 

Recall that for finite A^ T , a T and m r one may write in a toy model approximation plj] 

-5^ = ^ In det AlA Q = ^S;S Q , (40) 

a a 

where 

S Q (m') = cosh (AT T arcsinh (m r a T +ma T )) + —e^D, a + -e~^f2*. (41) 

Thereby at finite a T we may expect only approximate periodicity in A. That evidently differs from 
the case of m r — > 00 , which leads to ( |3S| ) . 

5 Zero modes 

Among the non-perturbative properties of quarks, those associated with the topology of the gauge 
fields occupy a special place. An important insight can be gained about the topological content 
of gauge field configurations in finite temperature QCD by computing the low lying modes of the 
Dirac operator. In the continuum the Dirac operator of massless fermions in a smooth background 
gauge field with non-zero topological charge has zero eigenvalues. The corresponding eigenfunc- 
tions are chiral and the number of left- and right-handed zero modes are related to the topological 



charge of the gauge configuration: ni — tir = v |2q] . Besides, the study of the Dirac operator 
spectrum near A = is of fundamental importance for the understanding of chiral symmetry 
breaking in gauge theories. 

Unfortunately, the lattice formulation also introduces artifacts and there is no exact index 
theorem in this context. Instead, what is observed with the actions currently used is the "zero 
mode shift" phenomenon |2^| (i.e. modes with small but non-zero eigenvalue and not exactly 
chiral). Typical gauge field configurations contributing to the path integral are not smooth and, 
in addition, lattice regularization breaks some other conditions of the index theorem as well: the 
topological charge of coarse gauge configurations is not properly defined and Dirac operator breaks 
chiral symmetry in an essential way pc| . Unlike the continuum, there is no division of lattice gauge 
fields into disconnected classes since lattice gauge fields are characterized by group elements on 
the links of the lattice and all link variables can be deformed to unity pi| . 

A noticeable progress in overcoming such difficulties was recently achieved . In particular, as 
it is argued in p5] |, the lattice QCD in fixed point formulation solves most of mentioned problems 
and in such context index theorem seems to be valid on the lattice. The lattice version of index 
theorem, which connects the small real eigenvalues of Wilson-Dirac operator and the geometrical 
definition of topological charge, was suggested in fl32| . 

Reliable zero mode separation in numerical computations still remains a serious problem []3lf , 
so we hope that even crude analytical models, which allow to estimate p(A) behavior in the area 
of small A, may appear useful. Notwithstanding the abrupt fall-off and vanishing of p(X) at A = 
measured in |33|] may be a finite volume effect |j) , the present data par excellence do not confront 
the suggestion that p(0) is very small and compatible with zero. As it was shown in (3^] the 
spectral density decreases as A approaches zero, and the larger is rif — the larger is the effect. 
The scenario when p(0) hits zero at nf > 5 is quite probable. 

Dirac spectrum should be particularly sensitive to the effects of quark loops since fermion deter- 
minant for small quark mass strongly suppresses gauge configurations with small Dirac eigenvalues. 
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Convincing arguments that spectral density disappears at A = are given in pM . As it is pointed 



out there, the weight in averaging involves a fermion determinant factor 

[Det(D - m)p = ' " 



A„>0 



(42) 



To be more specific we turn to the simple case of toy model approximation where fermion part 
of action can be totally separated from the gluon one and in accordance with ( |38|) incorporated 
into the measure by interchanging d/j,[ip a (x)] — ► dfi[ip a (x)} with 

N 

dp,[<p a (x)] = || (^cosh— + cos<y9 Q (x)J dfi[ip a (x)], (43) 

a=l 

with 

N / N \ N , 

n>m \k=l / k=l 

So the measure (H) has a set of n/-multiple zeroes, which evidently corresponds to the zeroes 
of (ft^). The effective action = SQ*[ip a (%)} for gluon field (e.g. the most popular Sq, 



suggested in |35|, 36 1) remains finite and scarcely able to influence either multiplicity or locus of 
the zeroes. Zero modes constitute no exception and are completely defined by the zeroes of (|4^). 
As it is easy to see from (p3), the eigenvalues of the matrix A x > yX , defined by det (A — iX) = 0, 



can be easily found from (^j) after the simple shift m — ► m — iX 

X = -im + T (ip a ± (2n + 1) n) (45) 
and one can easily show that the eigenvalues of the matrix A x , differs from ([i"5|) only by the sign 



of the real part. 

The expression ( |4q ) plainly shows that real modes appear only in the limit m — > and 
are placed at T (ip a ± (2n + 1) it) . It should also be noted that zero modes created by the fields 
^>(=fc) _ _^7o ^ are divideo; i n to pairs ±A as it is for = 1 ~ 75 ^> and j/)^ = 1 t i 75 -0. Although the 
total numbers of the pairs Ul + and n + + n_ are equal, the corresponding sets of eigenvectors, 
of course, do not obligatory coincide, so v — — hr and v' = n+ — n_ may differ. Therefore, 
the topological nature of j/, if any, needs special consideration. 
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It is not incurious that the spectral density can be treated simply as a result of integration of 
the QCD action over all variables <p x but one <p or the integration over the surface with the fixed 
average \ (<p) — y X (fx) ■ For example, in case of SU{2) gauge group one may write 

9 (f O) = nj (S (x V) ~ i £% &>-)) ) ^ = »/ exp {-S (1) ( V )} \ sin 2 (f ) . (46) 

The relation between A and the phases of Polyakov loop eigenvalues ip a makes it obvious that 
A are not fixed numbers, the dependence A on the parameters of the theory are defined not only 
by the fermion part of action, but by the gluon part as well. In particular, zero modes are defined 
by the relative number of Polyakov loops with cos ip a ~ — cosh S ~ — 1 . It is obvious that when 
the gauge fields are not fixed (<p a (x) ^ const) , the behavior of the corresponding det A (<p) in the 
limit m — > is generally different for different (p a (x) . Therefore, (fl6])does not result in an integer 
value for index v and we obtain the probability distribution, p(y), which is defined by (tp). 
However, we can compute the extreme value of v : i/ max = |max{z/}| . 

Taking into account the condition J2 a f a = '^ 7Tn we § e * "max = ^4 + ^j^j [^] , that both for 
SU(2) and SU(3) gives the same result ^ max = 4 + — . 

In a series of papers (37], |38|, |3^|, substantial new insight was added to this issue. Their 
central assertion is that the spectral density of Dirac operator very close to the origin A — should 
be universal, depending only on the symmetries in question. One of the amazing consequences of 
this conjecture is that the spectral density of Dirac operator near the origin need not be computed 
in full QCD, but can be extracted from simple models. 

We cannot exclude that simple analytical structure of Sp as in ([38]) will not persist in the 
general case, moreover, such simplicity is in conflict with numerical results obtained in JT^ ]. If, 
nevertheless, one ventures to apply universality arguments |37], |38|, [33] |3^] to our model, it will 
lead to p (A) ~ (^Y with v = v(g,T) and ^ max = 4+ Thus, the configurations with small 
eigenvalues are effectively suppressed, and the larger n/ is, the more prominent is the suppression. 
Therefore, in the computation of averages their contribution is small. It is true even for (i/'V') 
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because iptp has poles only in the points where fermion determinant turns to zero. 

6 Conclusions 

We have presented some evidence in favor of periodicity suggestion for the spectral density of 
Dirac operator on the lattice. Indeed, such periodicity appears in the limit a T — > in a toy model 
approximation |2l|| and in our opinion there are reasons to believe that it still persists in the 
general case. 

In this paper we take advantage of the fact that the chiral invariance of the massless part of 
Dirac operator leads to formal symmetry m <-> —m. This symmetry is of importance in studying 
QCD phase structure in (p; T; m) space ^0|. It appears desirable to make such symmetry explicit 
on finite size lattice. In case of infinitely heavy mirror fermion (r — 1), the symmetry m — m is 
restored. Moreover, the periodicity of spectral density has also become explicit for r = 1 already 
at finite a T . 

Besides, we suggest a simple way to relate the spectral density p (A) to the quasiaverage 
of QCD action with Polyakov loop matrices with arg f2 Q (x) = ip a (x) at some point x being fixed. 
As it can be anticipated (see e.g.,|y]]) the spectrum of Dirac operator is not invariant under Z(N) 
transformations and the results are different in the phases argfi = and argil = Nonetheless, 
for argfl = 4^ and argO = — the results coincide, which reflects C - symmetry (cp a <-> —<p a ) 
|p2| at zero density \i = 0. 

In this paper, we have taken a modest first step towards analytical investigation of Dirac 
spectral density characteristics. Apart from that there are many other questions which can be 
answered with the low lying eigenvalues and eigenvectors. The correlations between eigenvalues, 
the importance of which was proved in J3l| , is one of them. However, it is much more demanding 
a task than we can afford in such letter and we shall try to discuss it elsewhere 
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